Directional (version 3.1)

Density of the ESAG distribution: Density of the ESAG distribution

Description

Density of the ESAG distribution.

Usage

ESAGdensity(y, param, logden = FALSE)

Arguments

y

A matrix or a vector with the data expressed in Euclidean coordinates, i.e. unit vectors.

param

The parameters of the ESAG distribution, the first three are the mean vector in \(R^3\) and the next two are the two gammas.

logden

If you the logarithm of the density values set this to TRUE.

Value

A vector with the density values of y.

Details

MLE of the projected normal distribution, on the sphere, is implemented.

References

Mardia, K. V. and Jupp, P. E. (2000). Directional statistics. Chicester: John Wiley & Sons.

Paine P.J., Preston S.P., Tsagris M. and Wood A.T.A. (2017). An Elliptically Symmetric Angular Gaussian Distribution. Statistics and Computing (Accepted).

See Also

ESAGmle, ESAGsim, spml.mle, acg, circ.summary

Examples

# NOT RUN {
m <- colMeans( as.matrix( iris[,1:3] ) )
y <- ESAGsim(1000, c(m, 1,0.5) )
mod <- ESAGmle(y)
ESAGdensity( y, c(mod$mu, mod$gam) )
# }